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Abstract:  This  study  is  an  attempt  to  understand  the  flow  produced  by  a  hovering 
bird  or  by  a  hsh  when  it  starts  suddenly  in  stationary  water;  the  way  it  accelerates  and 
develops  thrust.  Major  focus  is  on  the  role  of  flexibility  of  the  flapping  foil  or  wing  in  these 
issues.  A  flexible  flap  with  negligible  mass  and  stiffness  is  attached  at  the  trailing  edge 
of  a  symmetric  airfoil.  The  airfoil  oscillates  sinusoidally  at  a  hxed  location  in  stationary 
water.  The  wake  flow  generated  by  two  airfoil  models,  one  with  flexible  flap  and  the  other 
with  rigid  trailing  edge,  are  compared;  the  former  produces  narrow,  coherent,  self-similar 
jet-like  flow  nearly  along  the  mean-position  of  model,  on  contrary  to  the  quite  divergent, 
meandering  flow  produced  by  rigid  trailing  edge  airfoil.  Comparison  suggests  that  the  flap 
plays  a  prominent  role  in  dehning  the  wake  signature.  Various  aspects  of  flap  motion  and  its 
flexibility  in  relation  with  the  flow  development  are  identihed.  Flexible  trailing  edge  induces 
multiple  vortices  along  with  large  vortices  arranged  in  ‘reverse  Karman  vortex  street.’  Flap 
accelerates  the  near-wake  flow.  It  sheds  the  vortices  with  appropriate  spacing  and  also 
imparts  them  convective  velocity.  Circulations  around  the  vortices  shed  from  flap  do  not 
increase  linearly  with  increase  in  the  chord.  Airfoil  with  flap  is  found  to  be  quite  robust  for 
amplitude-frequency  variation. 


1  Introduction 

The  pedantic  dehnition  of  ‘hovering’  suggests  that  the  mean  force  produced  by  the 
wings  is  exactly  vertical  and  precisely  balances  the  weight  of  the  creature  [8]  while  Rayner 
dehnes  hovering  as  the  mode  of  flight  in  which  the  body  is  at  rest  relative  to  undisturbed 
air,  with  the  wake  vertically  beneath  the  animal,  and  with  all  fluid  motions  induced  by  the 
beating  wings  [24].  To  achieve  this,  it  has  to  counter  the  gravity,  although  it  is  relieved  from 
the  act  of  generating  thrust  necessary  to  propel  forward.  The  animal  tunes  its  wing-flapping 
kinematics  such  that  thus  produced  force  will  be  just  enough  to  counterbalance  its  weight 
acting  vertically  downward,  enabling  the  bird  to  hold  the  desired  position.  To  fly  or  hover, 
momentum  needs  to  be  imparted  to  the  surrounding  fluid  which  they  do  in  the  form  of 
vortices  and  obtain  the  desired  forces  by  reaction.  Hovering  is  a  regular  activity  for  insects, 
but  for  birds,  it  is  “once-in-a-while”  affair.  So  insects  have  been  extensively  studied  to  delve 
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deeper  for  the  understanding  of  hovering  flight.  A  crucial  distinction  between  insect  and 
bird  flight  is  that  the  very  small  length  scales  involved  in  inset  flight  push  them  in  laminar 
flow  regime  whereas  almost  all  of  the  birds  fly  in  turbulent  regime,  generally  [11]. 

Three  types  of  hovering  are  identihed  depending  upon  the  angle  made  by  the  stroke- 
plane  with  the  horizontal.  The  first  type,  usually  known  as  ‘normal  hovering’  in  which  the 
stroke  plane  is  nearly  horizontal.  It  is  used  by  many  insects  and  also  by  the  hummingbird 
[11].  The  second  type  in  which  the  stroke  plane  is  inclined  and  makes  an  angle  between  30° 
to  60°  with  the  horizontal  is  known  as  ‘avian  hovering’;  e.g.  Dragonfly.  In  third  type  the 
stroke  plane  is  almost  vertical;  Ex.  Large  cabbage  white  butterfly.  The  lift  generation  takes 
place  either  by  circulation  or  prohle  drag  or  combination  of  both.  The  first  two  types  employ 
the  former  method  while  the  animals  hovering  in  the  third  type  rely  on  the  later  method 
[7].  For  the  later  two  methods,  the  force  is  mainly  generated  during  the  downstroke  while 
the  upstroke  produces  very  little  or  no  force.  Normal  hovering  mode  generates  nearly  equal 
force  in  both  the  strokes.  As  opposed  to  the  insects,  the  hummingbird  does  not  contribute 
equally  to  support  the  weight  during  both  halves  of  the  wing-beat  cycle  though  the  motion 
looks  kinematically  similar,  instead  75%  of  the  weight  snpport  comes  during  the  downstroke 
and  25%  during  upstroke  [31]. 

The  unusually  high  lift  coefficients  generated  dnring  hovering  are  quite  puzzling  and 
steady  as  well  as  quasi-steady  aerodynamics  theories  are  insufficient  to  explain  them.  Re¬ 
searchers  fonnd  that  hovering  animals  employ  unsteady  aerodynamic  mechanisms.  Wang 
[30]  and  Maxworthy  [21]  reviewed  the  nnsteady  aerodynamic  mechanisms  during  the  steady 
forward  flight,  gliding  and  hovering.  A  review  [25]  describes  the  fundamental  aerodynamic 
principles  lying  behind  the  insect  flight,  stressing  especially  on  the  unsteady  mechanisms 
involved  towards  the  refinement  of  the  existing  quasi-steady  model.  Sane  and  Dickinson 
improved  the  existing  qnasi-steady  model  by  adding  the  wing  rotation  effect  [26].  Amongst 
the  known  nnsteady  mechanisms,  some  are:  clap  and  fling  mechanisms,  delayed  stall,  vortex 
shedding  from  leading  edge,  wing  rotation,  rotational  circulation,  phase  adjustment  between 
wing  rotation  and  translation,  wake  capture,  flexibility  of  wings,  etc.  While  mimicking  the 
Dorsophila  that  hover  in  normal  mode,  with  robofly,  Dickinson  et.  al.  fonnd  that  the  high 
lift  coefficients  are  generated  by  three  distinct  and  yet  interactive  mechanisms:  delayed  stall 
dnring  translation  phase  while  rotational  circulation  and  wake  capture  generate  forces  dnr¬ 
ing  the  stroke  reversals.  The  phase  difference  between  wing  rotation  and  translation  is  also 
influential  in  force  generation  and  maneuvering  [6].  But  the  lift  enhancement  due  to  these 
unsteady  effects  is  accompanied  with  increase  in  drag  [30] .  Weis-Fogh  [32]  suggested  a  novel 
mechanism  of  lift  generation  known  as  ‘clap  and  fling’  for  hovering  insects  wherein  the  lift  is 
qnickly  generated  bypassing  the  Wagner  effect  [28].  Aerodynamic  analysis  by  Lighthill  [20] 
snggested  that  the  viscous  effects  and  the  hnite  span  of  the  wings  result  into  the  shedding 
of  vortex  rings  below  the  hovering  animal,  thus  supporting  the  weight  by  transforming  the 
momentnm  to  the  flnid  which  takes  the  form  of  a  jet  of  air  blowing  down  vertically  be- 


2 


low  the  animal.  Wang  [29]  obtained  nearly  similar  lift  for  a  two-dimensional  elliptic  wing  in 
hover  as  that  of  a  three-dimensional  wing,  suggesting  it  essentially  is  not  a  three-dimensional 
phenomenon.  Consideration  of  tip  vortex  modihed  the  ‘clap  and  fling’  model  [22]. 

To  smaller  or  larger  extent,  flexibility  is  an  ubiquitous  property  of  the  wings  and  hns 
of  flying  and  swimming  creatures.  The  wings  and  hns  often  deform,  quite  substantially 
sometimes,  when  they  propel  through  air  or  water  [5]  and  the  deformed  wing  or  £n  profile 
plays  key  role  in  determining  the  how  around  them,  making  the  analysis  even  more  profound. 
Due  to  deformation,  the  camber  of  the  wing  changes  continuously  and  slight  change  in  camber 
substantially  ahects  the  lift  generated  [1] .  Existing  literature  reveals  that  very  little  has  been 
contributed  towards  understanding  the  role  of  hexibility  often  found  in  the  force  generating 
appendages  which  are  partially  hexible,  except  a  few  studies  [14,  15].  Nachtigall  [23]  suggests 
that  among  other  mechanisms,  the  motion  of  hexible  wings  may  also  be  one  of  the  primary 
reason  for  the  large  lift  generation  during  hovering  [21].  The  ehect  of  wing  hexibility  during 
hovering  has  not  been  addressed  to  complete  details  due  to  the  intricate  huid-structure 
interaction  involved  [25].  That  makes  it  quite  occult  to  understand  analytically,  how  the 
wings  deform,  inertially  or  by  huid  dynamic  loading,  or  to  hnd  the  exact  contribution  of 
either.  Various  aspects  of  wing  hexibility  have  been  addressed  by  Daniel  and  Combes  [5]. 
The  aerodynamic  loading  is  quite  small  compared  to  the  inertial-elastic  forces  for  the  wings 
whereas  for  the  hns  in  water  the  huid  dynamic  loading  dominates  the  inertial  force  in  deciding 
the  hn  deformation  [5].  Wing  venation  in  insect  wings  determines  the  hexural  stihness  and 
its  spatial  distribution;  spanwise  and  chordwise  hexural  stihness  scale  respectively  with  cube 
of  the  wing  span  and  square  of  the  chord  length  [3,  4].  Zaho  et.  al.  observed  that  wing 
venation  along  with  wing  hexibility  is  very  important  in  force  generation  during  insect  hight 
|33|. 

Various  methods  have  been  used  to  analyze  the  hovering  hight  like  unsteady  actuator 
disk  theory,  jet-momentum  approach,  the  vortex  wake,  etc.  [24,  9,  10].  Ellington  used  vortex 
theory  and  modeled  hovering  with  “pulsed  actuator  disk”  that  applies  non-uniform  pressure 
impulses  to  the  air  [10].  The  happing  wing  of  hovering  animal  sheds  continuously  a  vortex 
sheet  in  sympathy  with  the  Kutta  condition  which  later  rolls-up  forming  the  vortex  rings. 
Rayner  analytically  modeled  the  wake  generated  from  rest  by  a  hovering  animal  with  a  chain 
of  coaxial  small-cored  circular  vortex  rings  stacked  one  upon  other  and  constructed  a  model 
to  predict  the  development  of  wake  shape  with  time  [24]. 

While  many  studies  are  available  for  thrust  generation  by  happing  foils  with  forward 
velocity,  very  few  studies  similar  to  hovering  have  been  conducted  wherein  the  happing  foils 
are  used  to  generate  thrusts  in  stationary  huid.  The  ability  of  a  sinusoidally  plunging  airfoil 
to  produce  thrust  is  known  as  the  “Knoller-Betz  ehect”.  Knoller  first  observed  this  fact 
in  1909  and  Betz  in  1912,  independently  [18,  2].  Katzmayr  observed  this  ehect  in  a  wind 
tunnel  in  1922  [17],  while  Jones  et.  al.  gave  its  experimental  and  numerical  proof  [16].  The 
dynamic  stall  vortices  can  be  fully  exploited  to  produce  large  thrust  coefficients  by  an  airfoil 
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in  hovering  and  the  hover  jets  can  easily  be  vectored  by  adjusting  the  initial  angle  made  by 
the  airfoil  chord  with  the  horizontal,  thus  allowing  easy  transition  from  hovering  to  other 
modes  of  flight  [12],  Lai  and  Platzer  [19]  observed  the  jet-like  velocity  prohles  downstream  of 
the  airfoil  plunging  at  zero  freestream  velocity  and  they  were  independent  of  the  amplitude 
and  frequency  of  oscillation  when  the  velocity  and  lateral  distance  were  non-dimensionalized 
by  the  maximum  plunge  velocity  and  amplitude  of  oscillation  respectively.  Heathcote  et. 
al.  [13,  14]  studied  the  effect  of  airfoil  flexibility;  the  moderate  spanwise  flexibility  increases 
thrust  as  well  as  efficiency  while  there  is  an  optimum  chordwise  flexibility  for  maximum 
thrust  generation  for  a  given  plunge  frequency  and  amplitude  combination.  The  flow  held 
was  independent  of  amplitude-frequency  variation  for  given  hexibility.  The  downstream 
momentum  hux  gave  the  rough  estimate  of  the  thrust  produced  by  the  airfoil. 

Another  issue  usually  encountered  in  the  thrust  generating  foils  is  the  jet  inclination, 
which  is  function  of  the  initial  conditions  and  random  jet-switching  which  could  be  triggered 
even  due  to  small  disturbances,  was  reported  earlier  [16].  Lai  and  Platzer  claimed  the  non¬ 
linear  interaction  of  closely  spaced  vortices  is  responsible  for  jet  inclination  [19].  Heathcote 
and  Gursul  [15]  showed  that  jet  switching  is  quasiperiodic  with  period  of  switching  two 
orders  higher  than  the  period  of  heaving  and  the  switching  frequency  varies  with  heaving 
frequency  and  amplitude  and  inversely  with  airfoil  stiffness.  Jet  switching  diminishes  with 
increasing  Strouhal  number,  irrespective  of  airfoil  stiffness. 

The  couplings  between  fluid  and  the  flexible  surfaces  are  not  yet  understood  completely 
[34],  so  we  approach  this  problem  experimentally  and  focus  to  answer:  How  important  is 
the  flexibility  often  observed  in  the  force  generating  flapping  elements,  in  these  issues?  This 
paper  investigates  the  effect  of  flexible  surface  appended  at  the  rigid  trailing  edge  on  the 
flow  over  a  symmetric  airfoil  pitching  sinusoidally  about  a  hinge  point  at  30%  chord  from  the 
leading  edge,  hxed  at  a  location  in  the  stationary  water,  in  an  attempt  to  mimic  the  normal 
mode  of  hovering,  usually  observed  in  insects  and  hummingbird.  The  flexible  flap  performs 
combined  heaving  and  pitching  and  traces  the  ‘figure-of-eight’  motion.  To  know  the  effect 
of  flexibility  and  the  difference  it  makes  to  the  flow  produced,  we  compare  it  with  the  flow 
over  another  airfoil  without  flap.  The  vortex  wakes  generated  by  the  two  airfoil  models  are 
studied  in  detail  in  a  plane  at  the  mid-span,  qualitatively  as  well  as  quantitatively.  The 
flows  produced  by  them  differ  considerably  conveying  that  the  role  of  the  flexible  flap  is  not 
insignihcant  even  though  its  mass  and  stiffness  are!  Flap  kinematics  is  studied  in  relation 
with  the  flow  generated.  Parametric  study  is  conducted  by  varying  amplitude  and  frequency 
of  oscillation. 

Some  of  the  probable  applications  where  this  study  may  prove  to  be  useful  are  outlined 
here.  The  flow  produced  during  the  current  experiments  is  nearly  similar  to  that  produced 
by  a  hand-fan  in  which  its  waving  creates  the  breeze  of  air.  This  may  be  useful  in  designing 
a  compact  device  that  produces  a  jet  of  fluid  by  simple  mechanical  oscillations  of  the  flexible 
flap,  towards  the  replacement  for  the  fan  in  electronic  cooling  applications.  Piezoelectric  fan 
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is  used  for  the  similar  purpose  where  the  leading  edge  vortices  are  used  to  create  the  flow. 
Though  in  this  paper,  we  are  not  dealing  with  the  analysis  of  unsteady  flow  produced  in 
the  very  beginning  i.e.  as  soon  as  the  flow  commences  during  the  current  experiments,  such 
study  will  provide  some  insight  in  understanding  the  force  generation  mechanism  employed 
by  the  fish  when  it  suddenly  accelerates  and  starts  to  move  from  rest  in  still  water-an  activity 
at  the  beginning  of  the  steady  swimming  and  also  during  unsteady  swimming,  common  in 
maneuvering  activities  like  catching  prays  and  escape  responses. 

2  Experiments 

The  experiments  are  conducted  in  stationary  water  in  a  glass  tank  (800  mm  X  800 
mm  X  350  mm).  The  airfoil-edges  in  spanwise  direction  are  placed  more  than  two  chords 
away  from  water  surface  and  bottom  wall  of  glass  tank.  The  model  is  confined  between  two 
end-plates  with  a  gap  of  about  1  mm  on  either  sides  as  shown  schematically  in  Fig.  1,  to 
keep  the  flow  two-dimensional.  Bottom  end-plate  is  made  of  glass  to  provide  optical  access; 
top  end-plate  is  of  Aluminium.  A  servo  motor  oscillates  the  airfoil  model  sinusoidally  about 
a  pitching  point  at  30%  chord  from  the  leading  edge.  The  desired  motion  is  prescribed  as, 

0  =  9maxSin{27lft),  (1) 

where,  ‘6*’  is  the  pitching  angle,  ‘F  is  time,  'Omax^  and  ‘/’  are  respectively  amplitude 
and  frequency  of  oscillation. 

It  is  a  Panasonic  A  series  AC  motor  (model:  RAMA3AZA1E),  30  W,  2500  ppr  in¬ 
cremental  encoder  controlled  by  a  closed  loop  servo  control  mechanism.  The  driver  is  a 
Panasonic  digital  AC  servo  drive  (model:  MSDA3A3A1A,  Type  1)  30  W,  220  V,  single 
phase  input  with  digital  and  analog  command  interfaces.  The  controller  card  is  DMC-1425, 
two  axis  card  from  Galil  Motion  Controls. 


Figure  1:  Schematic  of  the  experimental  set-up. 
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The  experiments  are  conducted  with  two  airfoil  models  (see  Fig.  2):  one  with  flexible 
trailing  edge  is  termed  as  model-F  and  the  other  with  rigid  trailing  edge  is  known  as  model- 
R.  We  fabricated  the  airfoil  using  the  ‘vacuum  casting  technique  The  airfoil  is  made  from  a 
plastic  material,  ABS  (Acrylonytril  Butadyne  Sulphide)  whose  density  is  1.01  gm/cc,  close 
to  water  density.  The  airfoil  section  is  NACA0015  symmetric  prohle  with  chord  of  the  rigid 
portion  40  mm  and  span  100  mm;  maximum  thickness  is  6  mm  at  30%  of  rigid  chord  from 
the  leading  edge.  We  extend  the  airfoil  chord  line  beyond  rigid  trailing  edge  by  appending 
a  flexible  flap  made  from  0.05  mm  thick  Polythene  sheet  which  has  negligible  mass  and 
stiffness.  Its  Young’s  modulus  (E)  is  3.02  x  10®  N/m^  and  flexural  rigidity  (El)  is  3.15  x 
10“^  N-m^.  It  is  inserted  in  a  0.1  mm  wide  slot,  that  protrudes  inside  the  airfoil  section 
for  10  mm,  and  provided  along  whole  span  at  sharp  trailing  edge  of  airfoil  (see  Fig.  2(c)). 
The  flap  exposed  to  the  fluid  is  30  mm  along  the  chord  and  100  mm  along  the  span.  Thus 
the  flexible  chord  is  75%  of  the  rigid  chord.  To  trace  the  flap,  it  is  blackened  in  the  plane 
of  visualization  except  near  the  trailing  edge  portion  to  identify  it.  Airfoil  is  provided  with 
two  dye  injection  holes  which  open  up  on  the  surface  at  the  maximum  thickness  location. 
Model-R,  the  airfoil  with  rigid,  sharp  trailing  edge  shown  in  Fig.  2(b)  is  exactly  identical  as 
model-F  except  that  there  is  no  flap  at  the  trailing  edge. 


Figure  2:  Schematic  sectional  views  of  the  airfoil  models:  (a)  model-F,  (b)  model-R,  (c) 
details  of  model-F.  All  dimensions  are  in  mm. 

The  flow  is  studied,  qualitatively  from  the  dye  and  particle  visualizations  and  quanti¬ 
tatively  from  the  particle  image  velocimetry  (PIV)  data,  in  a  horizontal  plane  located  in  the 
middle  along  the  airfoil  span  which  is  illuminated  by  a  laser  sheet.  The  optical  arrangement 
is  shown  schematically  in  Fig.  3.  The  laser  beam,  directed  through  the  front-coated  mirrors, 
is  converted  into  a  horizontal  sheet  using  a  combination  of  plano-concave  cylindrical  lenses 
along  with  a  plano-convex  spherical  lens.  We  captured  the  images  with  a  CCD  camera.  The 
distance  below  the  glass  tank  was  not  sufficient  to  focus  the  camera.  So  the  optical  path 
is  broken  and  the  visualizations  are  captured  through  a  flat  viewing  mirror  put  below  the 
tank.  Camera  unit  comprises  of  a  camera,  computer,  television-set  and  processor.  We  chose 
the  £eld-of-view  of  about  200  mm  x  150  mm  for  all  experiments  which  varies  to  some  extent 
from  experiment  to  experiment. 

For  flow  visualization,  we  use  a  continuous  laser  operating  at  up  to  4  W  at  488  nm 
(model;  Spectra-Physics,  Stabilite  2017).  The  visualizations  are  captured  through  Kodak 


(b)  (c) 


6 


Figure  3:  Schematic  illustrating  optics  arrangement. 

motion  analyzer  camera  and  Sony  camcorder.  We  conduct  dye  visualizations  with  fluorescein 
Sodium  salt  and  particle  visualizations  with  Polystyrene  particles  of  75-100  pm  diameter. 
To  track  the  flap  and  airfoil,  we  employ  shadow  technique  and  image  processing;  in  this,  the 
shadow  of  the  airfoil  and  flap  is  casted.  Airfoil  and  flap  locations  are  obtained  by  processing 
the  images.  The  error  introduced  during  clicking  and  picking  the  locations  is  about  1  to  2 
pixels  which  corresponds  to  about  0.15  to  0.30  mm.  Leading  edge  is  difficult  to  capture  for 
part  of  the  half  cycle  since  it  goes  into  the  shadow  of  the  airfoil  itself. 

The  velocity  measurements  are  carried  out  with  particle  image  velocimetry  (PIV).  In 
our  experiments,  the  near-wake  flow  is  two-dimensional,  so  the  out-of-plane  motion  is  negli¬ 
gible  there.  Hollow  glass  spheres  of  30  pm  size  have  been  used  as  tracer  particles.  We  use  a 
Nd:YAG  pulsed  laser  with  double  exposure  (model:  Big  Sky  Laser  Quantel,  ULTRA  CFR 
Nd:YAG  laser)  which  generates  120  mJ  of  energy  in  a  pulse-width  of  8  ns.  The  images  are 
acquired  using  Sharp  VISION  fire-wire  camera  (model;  1400-DE)  synchronized  with  pulsed 
laser  firing  which  is  controlled  through  a  computer.  All  images  spanning  the  fleld-of-view 
have  a  resolution  of  1360  x  1036  pixels  with  80  x  60  grid  nodes  in  X  and  Y  directions  respec¬ 
tively.  We  process  the  images  through  IDT  proVISION-XS  software  to  calculate  velocities. 
We  use  32  x  32  pixels  interrogation  window  size  with  nearly  50%  overlap  to  obtain  velocities 
at  all  grid  nodes.  Vorticity  is  calculated  using  a  subroutine  provided  in  the  software.  We  do 
the  data  post-processing  and  image  processing  in  MATLAB  6.0. 

The  parameters  varied  during  the  experiments  were  amplitude  and  frequency  of  flap¬ 
ping.  The  stall  angle  for  the  static  airfoil  case  is  about  14*^.  We  choose  three  amplitudes  of 
oscillation,  10°  at  which  there  is  no  static  stall,  15°  which  is  close  to  the  stall  angle  and  20° 
where  the  flow  may  separate  on  the  airfoil;  and  four  frequencies  for  each  amplitude,  1,  2,  3 
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and  4  Hz.  In  all,  we  collect  data  for  12  sets  of  parameters  listed  in  Table  1.  In  this  paper, 
we  present  a  detailed  analysis  for  Set-VI,  i.e.  amplitude  of  oscillation  15°  and  frequency  2 
Hz.  Except  otherwise  stated,  all  the  data  presented  in  this  paper  corresponds  to  Set-VI. 


Set 

I 

II 

III 

IV 

V 

VI 

VII 

VIII 

IX 

X 

XI 

XII 

amplitude 

(dynax  ) 

o 

o 

15° 

20° 

frequency 

(Hz) 

1 

2 

3 

4 

1 

2 

3 

4 

1 

2 

3 

4 

Table  1:  Values  of  the  parameters  used  in  the  experiments. 


3  The  mean  flow 

Sinusoidal  flapping  of  the  airfoil  with  the  flap  in  quiescent  fluid  creates  a  coherent, 
undulating  jet  composed  of  staggered  vortices  in  the  form  of  a  ‘reverse  Karman  vortex 
street’.  Figure  4  shows  an  ensemble  average  of  the  measured  velocity  and  vorticity  fields  at 
one  phase  during  the  oscillation  cycle  when  the  airfoil  trailing  edge  is  at  the  bottom-extreme 
location.  Although  the  jet  seen  in  the  mean  velocity  field  appears  to  move  nearly  along  the 
mean-position  in  the  downstream  direction,  the  ensemble  averaged  flow  shows  jet-meandering 
about  it  to  a  little  extent.  Dye  visualization  shows  (Fig.  5)  a  chain  of  smaller  vortices,  along 
with  two  relatively  large  counter-rotating  vortices  in  each  cycle.  These  smaller  vortices  are 
not  captured  in  the  PIV  data  due  to  insufficient  resolution.  Vorticity  plot  of  the  mean  flow 
held  in  Fig.  6  shows  a  different  view.  Two  strong  counter-rotating  vortices  draw  in  and 
accelerate  the  huid.  The  centers  of  these  vortices  are  near  the  hap  region.  The  shear  layers 
of  the  jet  portion  show  up  as  two  long  tongues  of  vorticity  of  either  sign. 
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Figure  4:  Ensemble  average  velocity  and  vorticity  field  when  the  rigid  trailing  edge  is  at 
bottom-extreme  location;  model-F. 
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Figure  5:  Dye  visualization  when  the  rigid  trailing  edge  is  at  bottom-extreme  location: 
model-F. 
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Figure  6:  Average  vorticity  from  mean  flow  field:  model-F. 

The  suction  of  the  quiescent  ambient  fluid  and  its  acceleration  to  form  a  coherent  jet 
are  clearly  seen  in  Fig.  7  which  shows  the  iso-contours  of  velocity  in  X-direction.  Such 
flow  behavior  is  explained  in  Sec.  4.  Jet-widening  due  to  entrainment  with  the  downstream 
distance  is  observed. 


-3  -2-1  0  1 


X/C 

Figure  7:  Iso-contours  of  velocity  in  X-direction  (u)  :  model-F. 

The  variation  in  magnitude  of  Umax  and  the  jet  half-width  with  the  downstream  dis¬ 
tance  are  also  plotted  in  Fig.  8.  The  entire  flow  domain  is  usefully  divided  into  different  zones 
in  X-direction  as  marked  over  the  figure.  The  maximum  u-velocity  shows  a  sudden  increase 
in  zone-u2  that  starts  from  near  the  mid-chord  location  of  airfoil  up  to  nearly  three-quarters 
chord  distance  downstream  of  the  rigid  trailing  edge;  this  is  the  active-flapping  region  where 
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work  is  done  and  momentum  is  added  to  the  fluid.  Umax  shows  steady  and  very  small  growth 
in  zone-u3  (from  -0.75C  to  -1.25C).  It  shows  a  rapid  increase  in  zone-u4  i.e.  beyond  -1.25C. 
The  highest  Umax,  which  is  136  mm/S,  is  attained  at  1.61  chords  downstream  of  the  rigid 
trailing  edge  and  close  to  the  mean-position  line.  Beyond  this  in  zone-u5,  a  steady  and  fairly 
gradual  fall  in  Umax  is  observed. 


Figure  8:  Variation  of  the  maximum  w-velocity  and  jet  half-width  in  the  downstream  direc¬ 
tion. 


Figure  9:  m- velocity  profiles  at  some  specific  downstream  locations;  model-F:  (a)  ^  ~  0,  at 
trailing  edge  (b)  ^  —0.76,  close  to  flap  tip  location  in  mean-position,  (c)  ^  ~  —1,  (d) 

^  —1.61,  highest  Umax  location,  (e)  ^  ~  —2,  (f)  ^  —3. 
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The  profiles  of  mean  u  (Fig.  9)  at  different  X-locations  show  near  Gaussian  variation 
except  at  the  edges.  The  fluid  flowing  with  high  u-velocity  is  confined  mostly  within  a  region 
spanning  half  chord  above  and  below  the  mean-position.  It  is  observed  from  Figures  8  and  9 
that  u  increases  in  downstream  direction,  even  beyond  the  flap  region,  for  certain  distance; 
this  increase  is  explained  by  the  actuator  disk  theory  in  Sec.  4. 

The  initial  narrowing  of  the  jet  width  is  associated  with  the  accelerating  flow.  Up 
to  two  and  half  chords,  the  flow  is  still  under  the  effect  of  the  favorable  pressure  gradient 
created  near  the  flap  region  (see  Sec.  4),  so  the  jet  grows  gradually.  Beyond  this  location, 
the  jet  grows  more  rapidly  due  to  entrainment  of  ambient  fluid.  Jet-width  is  nearly  one 
chord  at  downstream  distance  of  three  and  half  chords. 


4  Accelerating  flow:  actuator  disk  theory 

Figure  11  shows  the  streamlines,  and  pressure  and  velocity  variations  along  the  cen¬ 
terline  for  an  ideal  actuator  disk;  axial  momentum  flux  variation  will  be  same  as  that  of 
axial  velocity.  A  pressure  jump  across  the  disk  produces  favorable  pressure  gradients  and 
accelerating  flow  both  upstream  and  downstream  of  the  disk.  Momentum  increase  is  APA 
and  rate  of  work  done  by  the  disk  is  equal  to  (APA)|. 

The  variation  of  mean  momentum  in  X-direction  for  the  flapping  foil  experiment  (Fig. 
10)  is  similar  to  that  for  an  idealized  actuator  disk.  The  fluid  momentum  increases  slowly 
in  the  airfoil  region  (zones  ml  and  m2)  and  rapidly  in  zone-m3  where  the  flap  actively  flaps. 
It  then  increases  slowly  in  zone-m4,  rapidly  in  zone-m5  and  gradually  in  zone-m6.  momx 
remains  nearly  constant  thereafter,  but  with  very  small  fluctuations. 


Figure  10:  Variation  of  momx  in  downstream  direction  for  model-F. 

In  the  region  between  -1.75C  to  -2.5C,  Umax  keeps  reducing  (Fig.  8)  while  momx 
increases  (Fig.  10).  The  jet  widening  with  more  fluid  entrainment  and  effectiveness  of  the 
favorable  pressure  gradient  generated  near  the  model  in  increasing  the  velocity  of  entrained 


12 


fluid,  explain  this  flow  behavior. 


(b) 


Figure  11:  Schematic  of  the  ideal  actuator  disk:  (a)  the  flow,  (b)  velocity,  (c)  pressure. 
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Figure  12:  Streamlines  around  model-F,  from  the  mean  flow  held.  The  dotted  lines  above 
and  below  the  mean-position  line  indicate  the  jet-envelopes. 

The  accelerating  how  observed  in  the  immediate  downstream  (up  to  two  and  half 
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chords  downstream  from  rigid  trailing  edge)  of  model-F,  where  the  increase  in  momentnm  as 
well  as  velocity  (np  to  1.6  chords  downstream;  see.  Figures  8  and  9)  in  downstream  direction 
is  observed,  can  be  explained  by  the  actuator  disk  theory.  The  remaining  flow  domain,  where 
the  momentum  stays  nearly  constant,  can  be  modeled  like  a  jet. 

Streamlines  from  the  mean  flow  held  plotted  in  Fig.  12  show  that  the  mean  how 
generated  by  model-F  resembles  to  the  one  produced  by  an  actuator  disk.  There  are  some 
diherences  between  an  ideal  actuator  disk  and  model-F.  First,  in  the  idealized  case  velocity 
in  the  far  downstream  region  and  the  mass  hux  remain  constants.  In  contrast,  due  to 
entrainment,  mass  hux  increases  and  centerline  velocity  decreases  beyond  about  ^  =  —1.6; 
in  this  region  how  is  similar  to  that  in  a  two-dimensional  submerged  jet.  Till  about  ^  =  —2.5 
momentum  increases  due  to  the  favorable  pressure  gradient,  though  entrainment  is  present. 

The  other  diherence  come  from  the  fact  in  the  idealized  case  actuator  disk  is  an  in- 
hnitesimally  thin  disk  and  the  how  is  steady,  whereas  in  the  happing  foil  case  pressure 
increase  would  be  over  a  hnite  distance  and  the  how  is  periodic  in  time  due  to  the  sinusoidal 
pitching  of  the  foil.  Hence  the  pitching  airfoil  with  the  hexible  hap  may  be  more  suitably 
modeled  as  an  “unsteady  actuator  disk”,  in  which,  the  pressure  jump  and  the  velocity  in  the 
slipstream  are  the  functions  of  time.  The  pressure  diherence  which  draws  in  huid  and  pushes 
it  out  is  created  by  the  hap  motion.  When  the  model  oscillates,  the  hap  bends  considerably 
and  its  major  portion  keeps  pushing  the  huid  in  downstream  direction,  for  most  of  the  time 
over  a  cycle.  A  higher  pressure  would  exist  on  the  side  of  the  hap  which  is  pushing  the  huid 
whereas  its  other  side  would  experience  pressure  reduction  at  the  same  time.  Thus,  pushing 
action  of  the  hap  is  the  main  cause  for  the  speculated  sudden  pressure  rise  around  the  hap 
region. 


5  Conclusions 

We  studied,  qualitatively  and  quantitatively,  the  how  generated  by  two  airfoil  models, 
one  with  hexible  trailing  edge  i.e.  model-F  and  the  other  with  rigid  trailing  edge  i.e.  model- 
R,  when  oscillated  sinusoidally,  at  a  hxed  location,  in  stationary  water.  This  study  tries 
to  identify  the  role  of  hexibility,  often  observed  in  the  bird-wings  and  hsh-tails,  during  the 
how  generation  by  a  hovering  bird  or  by  a  hsh  when  it  starts  suddenly  in  stationary  water; 
the  way  it  accelerates  and  develops  thrust.  This  is  a  simple  method  of  how  generation  in 
stationary  huid,  in  the  sense  that  the  jet  and  hence  thrust  is  generated  just  by  the  mechanical 
oscillations  of  a  hexible  hap.  Although,  both  the  models  produce  jet-like  downstream  how 
sustained  mainly  by  the  vortices,  the  hows  diher  in  some  terms.  Addition  of  the  hexible 
hap  at  the  rigid  trailing  edge  of  airfoil  produces  more  streamlined  how  by  pulling-in  the 
huid  from  sides,  mainly  in  the  hap  region,  and  pushing  it  back  than  the  airfoil  without  it. 
Airfoil  with  hexible  hap  produces  a  narrow,  coherent,  undulatory  jet  conhned  among  the 
vortices  arranged  in  ‘reverse  Karman  vortex  street’,  whereas  the  one  without  hap  produces 
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quite  widespread  flow.  The  jet  produced  by  model-F  runs  nearly  along  the  mean-position 
with  slight  meandering;  on  the  other  hand,  the  jet  produced  by  model-R  meanders  about 
the  mean-position  quite  randomly  and  continuously.  We  attribute  this  jet  meandering  to 
the  vortex-vortex  interactions,  their  strengths,  shedding  pattern,  spacing  among  them,  etc. 
Interestingly,  almost  all  the  vorticity  is  conflned  within  the  narrow  region  around  the  mean- 
position  in  the  wake  of  model-F  whereas  it  is  spread  in  the  wake  of  model-R. 

Some  aspects  of  flap  kinematics,  studied  in  relation  with  the  observed  fluid  flow,  sug¬ 
gest  that  the  flap  alters  the  flow  substantially  and  plays  a  prominent  role  in  defining  wake 
signature.  Flap  tip  and  trailing  edge  move  almost  in-phase  in  X-direction,  but  a  phase  dif¬ 
ference  of  nearly  137°  exists  in  Y-direction;  trailing  edge  leads  the  flap  tip.  The  frequency  of 
flap  tip  deflection  in  X-direction  is  double  of  that  in  Y-direction.  Large  flap  deflections,  and 
even  second  mode  at  some  phases,  are  observed.  Flap  velocities  higher  than  rigid  trailing 
edge  affect  the  vortex  shedding  pattern.  While  model-R  mainly  sheds  large  vortical  struc¬ 
tures,  the  flexible  trailing  edge  sheds  vortices  of  two  kinds:  smaller  multiple  vortices,  shed  in 
order  to  satisfy  the  Kutta  condition  at  every  instant,  and  the  relatively  bigger  vortices.  The 
larger  vortex  is  the  manifestation  of  ‘stopping- vortex’  developed  during  the  direction-reversal 
of  rigid  airfoil  portion  which  is  strengthened  and  matured  while  moving  over  the  flexible  flap 
and  shed  by  its  tangential  removal.  The  rigid  airfoil  portion  plays  the  role  of  generating 
bigger  vortices,  whereas  the  flap,  apart  from  shedding  the  smaller  multiple  vortices  which 
may  later  merge  in  the  bigger  vortices  strengthening  them,  enhances  their  strength,  releases 
them  at  right  places  and  imparts  convective  velocity  due  to  which  they  and  hence  the  jet- 
like  flow  are  sustained  far  downstream.  Some  humps  seen  in  the  flap  tip  velocity  plot  in 
Y-direction  and  the  non-identical  flap  tip  acceleration  and  deceleration  in  X-direction  when 
it  travels  both-ways  might  be  the  consequences  of  shedding  bigger  vortex.  Moreover,  it  is 
shed  predominantly  only  when  the  flap  tip  moves  away  from  the  rigid  trailing  edge.  The 
jet-like  flow  is  found  to  be  stable,  qualitatively,  for  the  observed  vortex  aspect  ratio  (0.3584 
for  Set- VI). 

Flap  accelerates  the  flow  in  the  near-wake  region.  We  ascribe  this  to  the  favorable 
pressure  gradient  generated  due  to  the  sudden  pressure  difference  developed  by  the  flexi¬ 
bility  and  the  pushing  action  of  the  flap,  mainly  near  the  flap  region,  which  is  explained 
by  modeling  the  airfoil  with  flexible  flap  as  an  actuator  disk.  Momentum  in  downstream 
direction  continues  to  increase  even  after  Umax^  which  occurs  nearly  along  the  mean-position 
line,  started  to  reduce,  due  to  the  entrainment  causing  increased  fluid  mass  that  accelerates 
simultaneously.  Beyond  two  and  half  chords,  the  flow  obeys  characteristics  of  jet-like  flow. 
Total  chord  of  model-F  is  75%  more  than  model-R,  but  the  circulation  around  vortices  shed 
from  model-F  is  more  by  a  factor  of  about  three  than  model-R  which  is  a  consequence  of 
flexibility  introduced  by  the  flap.  A  point  of  zero  velocity  (POZV)  exists  over  the  flap  portion 
when  it  deflects  in  second  mode.  It  is  considered  as  a  movable  hinge  point  which  travels  over 
the  flap  region  resulting  in  differential  flap  portions,  pulling  and  pushing  the  fluid  about  this 
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hinge  point  and  thus  enhancing  the  circulation  shed,  finally.  The  downstream  fiow,  at  least 
within  the  jet-width  region,  is  found  to  be  self-similar  for  the  local  Umax  and  jet-width. 

The  airfoil  with  flexible  flap  is  found  to  be  quite  robust,  qualitatively,  for  the  amplitude- 
frequency  variation  since  it  produces  nearly  similar  fiow  for  all  the  parameter-sets  studied. 
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